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§1. Introduction 



Determining the phase structure of QCD in the chemical potential (yu)- temperature (T) 
plane is an important issue to understand properties of dynamics of QCD. Lattice gauge the- 
ory is one of promising tools to study non perturbative properties of QCD. At finite values 
of n, however, it suffers from notorious sign problem. Fermion determinant takes complex 
values and this invalidates importance sampling in Monte Carlo (MC) simulations. In order 
to circumvent the problem, various methods such as Taylor series method,'^''^ reweighting 
methocP^ etc. are attempted. Among such methods is analytic continuatiorP^''^''^''^''^''^ 
that we pay attention to in the present paper. In this method, one calculates some quantities 
in imaginary chemical potential region, fit the data by an appropriate function and extrapo- 
late it to the real /i region. Such an extrapolation, in general, requires careful assessment of 
extrapolated function. In order to explicitly check the validity of extrapolated functions, Cea 
et.al.liSKIU studied two-color qcDp^^^^^'^^^^'^^^^^^^^^^'^^ which is 
free from the sign problem. They calculated several quantities such as chiral condensate 
and Polyakov loop in the imaginary fi region, and compared these extrapolations with the 
results of direct computations in the real /i region. Based on what they studied, they made 
useful suggestions for applying the analytic continuation method to three-color QCD. In the 
current stage, however, it is a hard task to check their statements by using MC simulations 
because of the sign problem. It would then be useful to study these properties by adopt- 
ing some other methods. In the present paper we use chiral random matrix theory (RMT) 
modelfSSJ 

Chiral RMT is equivalent to QCD in the e regime. This microscopic equivalence is made 
use of for calculating low energy constants such as pion decay constant. From a mean field 
point of view, on the other hand, it is used to study phase diagram in the yU — T plane and 
has predicted relevant properties.'^^ Here, taking the latter stand point we study analytic 
continuations of two-color QCD in terms of the corresponding chiral RMT by investigating 
transitive regions from imaginary to real /i region at finite T. Based on this study, we shall 
investigate three-color QCD, the results of which will be reported in a forthcoming paper 
In the present paper, we adopt a model proposed by Klein, Toublan and Verbaarschot.'^''^ 
Since we are concerned about analytic continuation properties, we study the phase structure 
in the imaginary /i region, mainly in the vicinity of /i = 0. Although the phase structure 
in the real /i region is investigated in Ref. |32j), we present a part of these results associated 
with analytic continuation. 

It is found in MC studie^^^"^ that suitable fitting function depends on the region 
of temperatures, possibly reflecting the difference of involved physics. In particular, at 
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moderately high temperature region, the ratio type fitting function works well, and it is 
suggested to use it in the study of three-color QCD. In the present paper, we study the 
phase structure of the chiral RMT model in the real and imaginary n regions in the vicinity 
of /i = 0. Based on the results, we divide the temperature into three regions, T < (region 
I), To < T < Tco (region II) and T^o < T (region III), where is the maximal temperature 
of the diquark condensate phase in the real fi region, and Tco is temperature of cross over at 
^ = 0. We study the analytic continuation of the chiral condensate and the cross over line. 
As fitting functions, we use the polynomial types as well as the ratio type. In region III, we 
found that the polynomial fits show slow convergence, while the ratio type fits work well, 
in accordance with the results of MC simulation studies.'^ This is the main results of the 
present paper. 

The adopted RMT model incorporates a temperature effect only through the lowest 
Matsubara frequencies. This approximation is expected to be poor in the low temperature 
region. In order to check if the results stated above concerning the analytic continuation 
is sensitive to this approximation, we also consider an extended RMT model,'^ which 
incorporates all the Matsubara frequencies. We study the phase structure of this model in 
the real and imaginary /i regions. In the real n region, we make it clear the relationship of 
temperature with that in the original RMT model. In the region III, we found that the ratio 
type of function fits well the behaviors of the chiral condensate in this model, too. 

In the following section, we present the formulation of two-color QCD and its correspond- 
ing chiral RMT model. We briefly review its phase structure in the real fi region and study 
the phase structure in the imaginary fi region. We also present the results of the phase 
structure in the presence of a diquark source rather than quark mass. In section 3, we study 
analytic continuation of the chiral condensates and the pseudo critical line. In section 4, we 
discuss an effect of the higher Matsubara frequencies. Summary is presented in section 5. 

§2. Phase structure 

2.1. Two-color QCD and RMT 

In this subsection, let us briefly remind us of the chiral Lagrangian study and chiral 
random matrix model in two-color QCD for a completeness. Two-color QCD is characterized 
by pseudo reality 

T2raT2 = -T*, (24) 

where are Pauli matrices. As a consequence, the Lagrangian with massless Nf flavor 
quarks possesses enlarged Sl]{2Nf) flavor symmetry rather than SU(A^/)xSU(A^/) symmetry. 
An introduction of the chemical potential induces symmetry breaking SU(2A^/)^ SU(A^/)x 
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SU(A^/)x U(1)b. For non zero quark mass m 7^ 0, SU(2A^/) turns to Sp(2A^/), and further 
Sp(2A^/) SU(A^/)yX U(1)b due to a /i term. Some time ago, Kogut et. alP'^^ studied 
the model based on the chiral Lagrangian. The static terms to lowest order in the momentum 
expansion of the effective Lagrangian are given by 

C^^^^^ = -F^m^ReTr (^EM^ - F^^Tr {EB^E^B + BB) , (2-2) 

where is a 2Nf x 2Nf unitary antisymmetric matrix. The mass matrix M and baryon 
charge matrix B are given by 





where 1 is Nf x Nf unit matrix. The constant F is a low energy constant, and pion mass 
is given in terms of chiral condensate {ipip)o at = by = m{il)il))Q/{2NfF'^). The 
Lagrangian can be minimized by using the following form 

17 = l^o- cosa + 17^ sin a, (2-4) 

where 

For fi < mTr/2, the chiral condensate dominates, while for fi > m^/2 a rotation of condensates 
to the diquark condensate occurs as 

{Err < 'm-n/2) 

Z'o- cos tto + ^zi sin ao (/U > ""^77/2), 

where cosoq = m^/(4yU^). 

The QCD partition function with two flavors 

2 

^QCD = (nde*(^ + ^/ + mo)) (2-7) 
/=i 

is replaced by random matrix partition function 



JRMT 



VW e-x.Y>[--G'^TiW^W]detD, (2-8) 



where due to the pseudo reality, W is represented by a real x n matrix, and the probability 
distribution of the matrix elements is Gaussian. The matrix D is given bj^^ 

(mi iy + cu(T)+/ii \ 

ma W + u){T)+^i2 

-W'^-u{Tf + ^i mi 

-W^-uj{TY + H2 m2 I 
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As a temperature effect, only the Matsubara lowest frequencies are incorporated in the form 



T 
-T 



(2-9) 



In the present paper, we consider symmetric case fii = ^2, and mi = m2. 

Standard manipulations lead to a partition function represented by complex 2Nf x 2Nf 

{Nf = 2) matrix 



Zrmt = J VAexpl-C], 



C = -G^TiAA^ - ^ logdet Q+Q_, 



n 



(2-10) 
(2-11) 



where 



At + M fiB±iT 
-fiB ±iT A + M^ 
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(2-12) 



and B = diag(l, —1, 1, —1). 
An ansatz for A 
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(2-13) 



leads to an effective potential: 
1 



n 



C = 2G\a' + A')-\yiD, 



(2-14) 



D= [[a + m + ^iY + A'^ + T^] [{a + m - ^if + A'' + T^] . 

2.2. saddle point equation 

In order to obtain the most dominant contribution in Eq. fl2-10p in the n — 00 limit, we 
solve the saddle point equations given by 

dC 



This yields 



aG^D - ((a + mf - fi^ + A^ + T^) (a + m) = 0, 
A [G^D - {{a + mf + + A"" + T^)] = 0, 



(2-15) 

(2-16) 
(2-17) 
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where D reads 



m] 



(2-18) 



Saddle point equations in the imaginary n region are given by replacing fi by in Eq's. fl2-16p . 
(12 -171) and fl2-18p . In the following subsections, the phase structure in the imaginary /i region 
is studied. In order to understand the following analysis of the analytic continuation, we also 
discuss the phase in the real /i region. Although the latter is studied in detail in Ref. 1521) . 
we present a part of it in a self-contained manner. It should be noted that the RMT model 
does not show a periodicity in the imaginary /i, caused by the center symmetry of the gauge 
group (Roberge- Weiss symmetrj^® ). This point is discussed in summary in connection 
with the analytic continuation. 

2.2.1. m = 

In the following section, we shall discuss the analytic continuation in the case of non zero 
quark mass. In order to have a better understanding of the phase structure in the massive 
case, we also discuss the phase structure in the massless case. For m = 0, the saddle point 
equations are analytically solved. For T = 0, the saddle point equations become 

a = 0, G^D = a^~fi^ + A^, , , 

f2-19) 

Z\ = 0, G^/} = fT^ + + ^ ^ 



D = ((T^ - + Z\^)^ + 4Z\V- (2-20) 

This is classified into the following four cases:(a) a = 0, Z\ = 0, (b) a = 0, Z\ 7^ 0, (c) 
0" 7^ 0, Z\ = 0, (d) 0" 7^ 0, Z\ 7^ 0. Correspondingly, the effective potential reads 

Q = -C^"^ = -21n/i2, 

n 

^(A) ^ 1^(6) ^ 2(1 + In G2 _ ^2^2) (^2 ^ ^2 ^ ^/^2)^ 

n 

Qi^) ^ l£(c) = 2(1 + InG^ + GV') - /i' = 
n 

= i/:^ = 2(l + lnG'2) at/i = Oonly. (2-21) 

n 

For < /i < /Xc, i^^"^^ gives the lowest values among these i7, i.e., diquak condensation 
occurs, where /ic is determined so that 2(1 — G'^fj?) + ln(G^/i^)^ = (i?^'^) = Q) is satisfied. 
For /i > /ic, symmetric phase realizes. For imaginary = icf), one can easily see in Eq. fl2-2ip 
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that minimum is found just by an interchange of a and Z\, i.e., the roles of cr and A are 
interchanged in the real and imaginary /i regions. For < < 0c, chiral condensation 
occurs, and for 0c < 0, o" = 0, Z\ = 0. The value of 0c satisfies 2(1 - C^cfl) + \n.{G'^ (f^Y = 
= where i?*^'^-' = H^'^^fi — > i0) etc.. At = 0, the system undergoes a first 

order phase transition. When one crosses the origin from real /i to imaginary /i region, the 
magnitudes of A in the real /i region jumps down to at = 0, while that of a jumps up 
to a finite value from vanishing one. 

At finite temperatures, this first order critical point persists up to a finite value Tc along 
the T-axis as shown in the left panel in Fig{Tl Temperature effect also shifts the location of 
/ic such that 



l + \nG\fil + T^)-G\fil + T^ 



0, 



while 0c is determined by 



2G\<j)l + T^) - (l + 



1602r2G4j = 0. 



(2-22) 



(2-23) 



The second order critical line in the real fi region starts at /i = fic, passes the endpoint 

2n 



T 
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■I -1.5 -1 -0.5 0.5 



Fig. 1. Phase diagram for m = (left) and m = 0.1 (right). Solid lines indicate second order 
critical lines. Dashed dotted and dotted lines indicate first order critical line and cross over, 
respectively. G = 1.0. Note that negative side of the horizontal axis indicates the imaginary 
chemical potential, 



{fi = 0, Tc) of the first order critical line at /x = and extends to the imaginary region as T 
increases. On the other hand, the critical line, separating symmetric and chiral condensate 
phases in the imaginary /i region, starts at 0cG = 1.0 and T = 0, and moves toward larger 
values of as T increases. These are shown in the left panel in FigUl 
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2.2.2. m^O 

For m 7^ 0, the effect of chiral condensate becomes stronger and the phase structure 
changes. The first order critical hne, separating chiral and diquark condensate phases at 
/z = in the m = case, changes to the second order one and shifts towards real /i region as 
shown on the right panel of Fig{TJ This is given as the lowest energy solution of Eq's. fl2-16p 
and fl2-17l) for small values of fi. At T = 0, the phase transition occurs at ficG = 0.26 = ficiG. 
For fi > fici, chiral condensate rotates into diquark condensate, which is in agreement with 
the chiral Lagrangian study as in Eq. (12-61) . Figure [2] indicates contour plots of the effective 
potential in the a — A plane for m = 0.1 and T = {fiG = 0.0,0.4 and 0.8 from left to 
right). The darker the plot becomes, the smaller the magnitude of the effective potential 
becomes. We see that the minimum rotates from [a 0, A = 0) to {cr = 0, A ^ 0) as fi 
increases. Unlike the chiral Lagrangian, however, diquark condensate becomes vanishing at 
fiC = ficG{= 1.0) at T = 0, and the diquark condensate phase turns to quasi symmetric 
one (cr ^ 0{m),A = 0). This is due to the saturation effect coming from the limited 
dirac spectrum of the RMT. As temperature increases, diquark condensate phase becomes 
narrower and terminates at (/x^i. To). For T > To, thermal fluctuations become large enough 
to make diquark condensate cease. A cross over line leaves the point {^DiTd) and moves 
towards /i = as temperature increases, where the signal of cross over is detected by a peak 
of chiral susceptibility. This line is a reminiscent of the second order critical line in the m = 
case. In the imaginary region, it further proceeds toward larger values of (f) as temperature 
increases. For larger values of 0, the signal of the critical line in the m = case, emanating 




-2-10 1 2 -2-10 1 2-2-1 1 2 

Fig. 2. Contour plots of the effective potential va. a — A plane for m = 0.1 at T = 0. From left to 
right, ijl = 0.0,0.4 and 0.8. The cross symbol indicates the location of the minima. G = 1.0. 
The horizontal (vertical) axis stands for a (Z\). 

at 0c, becomes so weak in the m 7^ case that the peak of the susceptibility is very broad 
and low. In the figure, we do not show it. 
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2.3. diquak source jd 7^ 

So far we have discussed the phase structure in the case of non zero quark mass. In 
two-color QCD, meson and diquark states belong to the same representation, and then 
mass and diquark source are on an equal footing. It is interesting, therefore, to discuss the 
phase structure when a diquak source is switched on. When m = 0, a diquark source jd is 
introduced in Eq. f l2-12p as 



M 



The effective potential becomes 

1 
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/ 



(2-24) 



-C 

n 



\2 



1 2 



Jd) 



[{A + jdf + a' - ^^'Y + 4f^\A 
+ 2[{A + jdf + a^ + /x^] + T\ 

At T = 0, the argument of the logarithmic term in Eq. (12-251) reads 

D^'^=[{A + Jdf + a^-^?]' + A^^\A+Jdf 



(2-25) 
(2-26) 

(2-27) 



Comparing this with D in Eq. fl2-18p at T = 0, we recognize that there is a symmetry with 
respect to a simultaneous interchange of jd ^ A ^ a and /i^ ^ ■ This implies that 
the behavior of a {A) in the 7^ 0, m = case is a mirror image of the behavior of A (cr) in 
the jrf = 0, m 7^ case if one plots them as functions of /i^ as shown in Fig. [3l The rotation 
of condensate thus occurs in the imaginary /i region, from Z^7^Otoo"7^Oas0 increases 
from 0. 

At finite T, a term proportional to breaks this symmetry as seen in Eq.'s (12-18P and 
fl2-26l) . At low temperatures, such a breaking effect is small, but at higher temperatures, it 
becomes clearer. By investigating the behaviors of a and A in the similar way to the m 7^ 
case, we obtain the critical lines as shown in Fig. HI Unlike the m 7^ 0, = case, the phase 
with both of a and A condensation, located now in the imaginary region, does not close 
itself but stretches to larger temperature region. 

§3. Analytic continuation 

In this section, we discuss analytic continuation for m 7^ 0. Hereafter, G is taken to be 
1.0. Figure [5] indicates free energy and chiral condensate in the real and imaginary /i regions 
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1 .5 



Fig. 3. a for (j,m) = (0.1,0) and A for (j,m) = (0,0.1) at T = 0. G = 1.0. Negative side of the 
horizontal axis indicates the imaginary chemical potential, 
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Fig. 4. Critical lines for j = 0.1 and m = 0. G = 1.0. Solid line indicates the second order critical 
line. Negative side of the horizontal axis indicates the imaginary chemical potential, 



for various temperatures for m = 0.1. Both of the quantities behave smoothly at /i = for 
all temperatures considered here. Based on these results, we check analytic continuation in 
terms of the RMT and compare the results with those of MC simulations.'^^"^ 

3.1. chiral condensate 

We use calculated values of a in the imaginary fi region as "data" . We take s points of 
data in a certain range < < 0/, fit them by an appropriate function and extrapolate it 
to the real region. Following Cea et. al., the following two types of fitting function are 
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-0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 -0.5 
.2 



Fig. 5. Free energy (left) and chiral condensate (right) in the vicinity of fi = for m = 0.1 and 
various values of T. 



used, one being polynomial type and the other ratio type. 



A + B(jP + C(j)'^ + ... (polynomial), 

(ratio). 



A+B(f>- 
l+C<t>' 



For a systematic study, we check to what extent the outcome depends on the following 
points: 

1. the range of used "data" in the imaginary /i region 

2. the number of degrees of freedom of the data 

3. the highest degree of the polynomial in the case of polynomial fit 

We used 0/ = 0.3, 0.5 and 0.8, and found that the outcome depends rather strongly on the 
adopted range 0/. Generally speaking, the 0/ = 0.3 case turns out the best among all the 
cases we investigated. On the other hand, when we varied the number of degrees of freedom, 
d.o.f., from 1 to about 75, depending on data, the results depend hardly on the number of 
data s. This is shown in Table H] for the T = 0.8 and 1.5 cases. Hereafter, we present the 
results of 0/ = 0.3. 

We divide temperature into three different regions, T < To (region I), To < T < Tco 
(region II) and T^o < T (region III), where is the maximal temperature of the diquark 
condensate phase (T^ = 0.892) and Tco is temperature of cross over at /i = {Tco = 1-18). It 
should be noted that this classification of the region of temperature are different from that 
in Ref . [TUl) . Regime a and b in Ref. | TU1 ) correspond to our region III, and regime c in Ref. [TU|) 
are further classified into region I and II in our case. In the latter, we distinguish whether 
analytic continuation hits the diquark condensate phase (region I) or not (region II). In region 
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Table I. Dependence of fitting parameters on the number of data for T = 0.8 and 1.5. The adopted 
range for the data is 0.0 < (p < 0.3. A(j) is defined by the number of points s as 4>f/s. 



f{<j)) = A + B<f)^ + C(j)^ 





T = 


0.8 


T = 


1.5 


Acj) 


0.01 


0.05 


0.01 


0.05 


d.o.f. 


28 


4 


28 


4 


xVd.o.f. 


5.63814 X 10-9 


1.08933 X 10-8 


5.2598 X 10-11 


1.02758 X 10-1° 


A 


0.62178(2) 


0.62178(7) 


0.07803(0) 


0.07803(0) 


B 


-0.68756(169) 


-0.68116(454) 


-0.17377(16) 


-0.17360(44) 


C 


-1.67541(2050) 


-1.6443(506) 


0.32745(198) 


0.33086(492) 


g{(j)) = A + B(I)'^ + Ccj)^ + D(p^ 




T = 


0.8 


T = 


1.5 


Acj) 


0.01 


0.05 


0.01 


0.05 


d.o.f. 


27 


3 


27 


3 


xVd.o.f. 


3.18334 X 10-11 


5.12716 X 10-11 


3.03625 X 10-1^ 


4.60497 X 10-1^ 


A 


0.62171(0) 


0.62170(0) 


0.07802(0) 


0.07802(0) 


B 


-0.70553(29) 


-0.70565(74) 


-0.17551(0) 


-0.17550(2) 


C 


-2.26163(849) 


-2.26256(2152) 


0.27069(2) 


0.27073(65) 


D 


-4.65104(6623) 


-4.63419(15920) 


-0.45033(205) 


-0.45081(477) 






M0) = ^ 


2 
2 






T = 


0.8 


T = 


1.5 


Acj) 


0.01 


0.05 


0.01 


0.05 


d.o.f. 


28 


4 


28 


4 


xVd.o.f. 


4.40095 X 10^11 


9.37848 X 10-11 


1.93257 X 10-12 


3.59884 X 10-12 


A 


0.62169(0) 


0.62169(0) 


0.07802(0) 


0.07802(0) 


B 


-2.88150(288) 


-2.88685(769) 


-0.05713(15) 


-0.05738(37) 


C 


-3.49352(434) 


-3.50174(1153) 


1.52002(163) 


1.51722(399) 



I, we do not consider very low temperature because this model incorporates only the lowest 
Matsubara frequencies and is thus expected to be poor in predictability. The left panel in 
Fig. indicates the results of a for T = 0.8. At this temperature, there exists a critical point 
at fill — 0.089. Extrapolations to the real fi region for both of the two types of fitting are 
good up to this critical point /id, but becomes extremely bad beyond /Xci, as it is expected. 
This behavior corresponds to the behavior in the low temperature region in Ref . [TOj) . When 
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they use a polynomial function as an interpolation at this temperature region, the deviation 
between the fit and the real data occurs at the point where a chiral condensate susceptibility 
peak occurs in the real jj region.'^ In our case, this peak is observed as a sharp one at 
yUci, where the diquark condensate starts to rise. In the polynomial case, if the higher power 
contributions are included, the above stated feature does not change (see Table [IT|) . 
In region II, we choose T = 1.0. At this temperature, a cross over exists at /i^ = 0.109. In 

Table II. Values of data and the results of various fits at real values of /x for T = 0.8. The range 
of used data is 0.0 < (j) < 0.3, and Acp = 0.01. Note that a critical point exists at /x ~ 0.299. 





fi = 0.0 


fi = 0.1 


fi = 


0.2 


fi = 0.3 


fi = 0.4 


/i = 0.5 


data 


0.63270 


0.61437 


0.5^ 


^886 


0.52881 


0.23333 


0.10833 


4th degree 


0.62178 


0.61474 


0.59160 


0.54633 


0.46888 


0.34518 


6th degree 


0.62171 


0.61442 


0.5^ 


^957 


0.53650 


0.43187 


0.23130 


8th degree 


0.62170 


0.61437 


0.5J 


^905 


0.53233 


0.40771 


0.12378 


10th degree 


0.62170 


0.61437 


0.5J 


^891 


0.53046 


0.39102 


0.01637 


ratio 


0.62169 


0.61434 


0.5^ 


^870 


0.52854 


0.36526 


-0.77936 



the polynomial method, the effect of the higher power terms becomes smaller, and the 4th 
degree of polynomial yields as good results as the higher degree of polynomials as shown in 
the right panel in Fig. [61 Deviations from the data become large at a value of fi slightly 
higher than the cross over point where a chiral condensate susceptibility peak occurs. The 
ratio method is also in the same quality as the polynomial ones. In Table IIIH we show the 
numerical values of the results. 

In the higher temperature region, region III, however, we obtained different results. Figure [7] 

Table III. Values of data and the results of various fits at real values of // for T = 1.0. The range 
of used data is 0.0 < (p < 0.3, and Z\(/> = 0.01. Note that a cross over point exists at ^ ~ 0.33. 





fi = 0.0 


= 


0.1 


fi = 0.2 


fi = 0.3 


/i = 0.4 


/i = 0.5 


data 


0.40000 


0.3^ 


^615 


0.34111 


0.25935 


0.16310 


0.09460 


4th degree 


0.40007 


0.34220 


0.38640 


0.25780 


0.11710 


-0.10244 


6th degree 


0.40000 


0.3^ 


^611 


0.34035 


0.24884 


0.08340 


-0.20614 


8th degree 


0.40000 


0.3^ 


^614 


0.34060 


0.25083 


0.09493 


-0.15482 


10th degree 


0.40000 


0.3^ 


^615 


0.34092 


0.25511 


0.13320 


-0.09155 


ratio 


0.39996 


0.3^ 


^593 


0.33880 


0.23825 


0.01903 


-0.62254 



indicates the fitting at T = 1.5. In the left panel, we compare the different types of fitting. 
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Fig. 6. Fitting at T = 0.8 (left) and T = 1.0 (right) for m = 0.1. Two types of polynomial fit (the 
4th and 6th degree) and ratio type are compared with the "data". 

and find that the ratio type of fit works very well. This behavior is in good agreement with 
that in the higher temperature region (regime a and b) in the MC simulations.'^ In contrast 
to this behavior, the polynomial fit shows the behavior with very slow convergence as shown 
in the right panel. There, we compare the results of the polynomial fit with the maximal 
power from 4 to 10. In Table HVl we show the numerical values of these fittings at several 
values of fi. In order to check if the ratio method works well at other temperatures in region 
III, we chose T = 1.7, 1.8, 1.9,2.0 and 2.2 as temperatures higher than 1.5, and found that 
this is indeed the case. When we chose, on the other hand, the value of T slightly lower 
than 1.5 in region III, the same tendency is seen, but the contrast between the ratio and 
polynomial types becomes less clearer than the case of T=1.5. The results at T = 1.3 are 
listed in Table |Vl The values of the fitting parameters used for T = 0.8, 1.0, 1.3 and 1.5 are 
summarized in Table IVTl 

3.2. chiral susceptibility peak 

The analytic continuation method has been applied to the determination of the critical 
line.Sl' We here focus on the chiral susceptibility peak for this purpose. The chiral suscepti- 
bility Xch is calculated by the m-dependence of a. The location of a peak of Xch in /i^ — T 
plane for m = 0.1 is shown in Fig. [8l We fit the imaginary "data" for < < 0.32 by 
polynomial formulae, one is A + B^^ and the other A + + C/x^ + Dfi^, and extrapolate 
them to the real region. These are compared with the results for real n values of the model 
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Table IV. Values of data and the results of various fits at real values of fj, for T = 1.5. The range 
of used data is 0.0 < (j) < 0.3, and Acf) = 0.01. 





fi = 


0.0 


/i = 0.2 


/i = 0.4 


/i = 0.6 


fi = 0.8 


/i = 1.0 


data 


0.07; 


m 


0.07142 


0.05570 


0.03837 


0.02399 


0.01353 


4th degree 


0.07; 


m 


0.07160 


0.05861 


0.05791 


0.10094 


0.23171 


6th degree 


omi 


m 


0.07141 


0.05503 


0.02891 


-0.04148 


-0.2771 


8th degree 


0.07; 


m 


0.07142 


0.05578 


0.04097 


0.05876 


0.26905 


10th degree 


0.07; 


m 


0.07142 


0.05573 


0.04058 


0.05345 


0.22534 


ratio 


0.07; 


m 


0.07139 


0.05540 


0.03713 


0.02101 


0.00829 



Table V. Values of data and the results of various fits at real values of fi for T = 1.3. The range 
of used data is 0.0 < cj) < 0.3, and Acj) = 0.01. 





/i = 0.0 


fi = 0.2 


fi = 0.4 


fi = 0.6 


fi = 0.8 


fi = 1.0 


data 


0.13425 


0.11678 


0.07971 


0.04646 


0.02424 


0.01072 


4th degree 


0.13426 


0.11721 


0.08920 


0.11963 


0.32415 


0.86470 


6th degree 


0.13425 


0.11694 


0.08420 


0.07916 


0.12544 


0.15474 


8th degree 


0.13425 


0.11671 


0.07362 


-0.09124 


-1.29117 


-7.56647 


10th degree 


0.13425 


0.11678 


0.08174 


0.16905 


2.37356 


22.91351 


ratio 


0.13422 


0.11625 


0.07496 


0.03110 


-0.00495 


-0.03183 




Fig. 7. Fitting at T = 1.5 for m = 0.1. Two types of polynomial fit (the 4th and 6th degree) and 
ratio type are compared with the "data" in the left panel. The maximal power dependence in 
the polynomial fit is shown in the right panel. 
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Table VI. Parameters of various fits of the chiral condensate in accordance to ct((/>) = A + Bcj)^ + 
C(t)^ + D(t)^ + E<f + F<i)^^ (polynomial; p) or ct((/>) = , (ratio; r). Note that blank columns 

stand for terms not included in the fit. 



T 


Fit 


A 


B 


C 


D 


E 


F 


d.o.f. 


xVd.o.f. 






P 


0.62178 


-0.68756 


-1.67541 








28 


5.63814 X 10 


-9 




P 


0.62171 


-0.70553 


-2.26163 


-4.65104 






27 


3.18334 X 10- 


-11 


0.8 


P 


0.62170 


-0.70780 


-2.39816 


-7.22101 


-14.9146 




26 


1.97889 X 10- 


-13 




P 


0.62170 


-0.70807 


-2.42373 


-8.05784 


-25.9669 


-50.6777 


25 


1.57648 X 10- 


-15 




r 


0.62169 


-2.88150 


-3.49352 








28 


4.40095 X 10- 


-11 




P 


0.40007 


-1.33932 


-2.68272 








28 


4.65810 X 10 


-9 




P 


0.40000 


-1.35568 


-3.21656 


-4.23547 






27 


8.22244 X 10- 


-12 


1.0 


P 


0.40000 


-1.35459 


-3.1514 


-3.00882 


7.11877 




26 


1.05268 X 10- 


-12 




P 


0.40000 


-1.35399 


-3.09274 


-1.08930 


32.4704 


116.244 


25 


2.02417 X 10- 


-14 




r 


0.39996 


-2.43151 


-2.66407 








28 


1.26412 X 10 


-9 




P 


0.13426 


-0.47440 


1.20483 








28 


1.08238 X 10" 


-10 




P 


0.13425 


-0.47683 


1.12564 


-0.62833 






27 


6.11875 X 10- 


-12 


1.3 


P 


0.13425 


-0.47782 


1.06582 


-1.75429 


-6.53443 




26 


4.66163 X 10- 


-14 




P 


0.13425 


-0.47769 


1.07827 


-1.35679 


-1.15240 


24.6779 


25 


5.27750 X 10- 


-17 




r 


0.134221 


-0.22694 


1.91334 








28 


9.61851 X 10- 


-10 




P 


0.07809 


-0.17377 


0.32745 








28 


5.2598 X 10- 


11 




P 


0.07802 


-0.17551 


0.27069 


-0.45033 






27 


3.03625 X 10- 


-14 


1.5 


P 


0.07802 


-0.17544 


0.27492 


-0.37034 


0.46189 




26 


1.51937 X 10- 


-17 




P 


0.07802 


-0.17544 


-0.27490 


-0.37132 


0.45424 


-0.03506 


25 


1.57037 X 10- 


-17 




r 


0.07802 


-0.05710 


1.52002 








28 


1.93257 X 10- 


-12 



as shown in Fig. [HI The latter fits very well, while the former does not function at all. That 
higher order terms play a relevant role in the case of chiral susceptibility peak is in good 
agreement with that of Cea et.al. (see Fig's. 3 and 6 in the first reference in Ref. ITTj) ). 

§4. Effect of higher Matsubara frequencies 

The model we have studied so far, let us tentatively call it the KTV model, incorporates 
temperature effects only from the lowest Matsubara frequencies. This limitation is expected 
to give an influence on the behaviors of the baryon density and chiral condensate at low 
temperatures, because the Matsubara frequencies form a continuum at zero temperature 
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1.4 



T 




0.15 



Fig. 8. Locations of the chiral susceptibility peaks in fi^ — T plane for m = 0.1, and two types 
of polynomial fit: linear in ijp, 1.180 — 1.171fj? (solid line) and the maximal power with fi^, 
1.178 - 1.294//2 _ 2.157//^ - 11.564/^^ (broken line). 

limit. Vanderheyden and JacksorP^ introduced an extended RMT model (VJ model) in- 
corporating contributions from higher Matsubara frequencies. We study this model in the 
imaginary n region and investigate the analytic continuation in order to check to what extent 
the results obtained in the previous section are robust, particularly, in region III. 
@ The partition function of the VJ model is given by 



dadAe 



-riQ 



(4-1) 



The effective potential is given by 

fI = 2l3G\a^ + A^)-f{a, Z\), 

oo 

/(a, Z\) = 5^ 5^ In ((a + m ± ^,f + Z^^ + (2n + ifT^'T^)] , (4-2) 

where the summation over n stands for the one over Matsubara frequencies, and the sum- 
mation over ± does their positive and negative parts, and (3 = 1/T. A term, referred to as 
the regular term in Ref. [33]), is dropped in Eq. ( ]4-2p . since it is independent of a and Z\, and 
irrelevant in the following argument concerning the analytic continuation. By defining 



the partition function reads 



dadAe 



-nl3Q 



(4-3) 
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where f2 is expressed as 

oo 

-P-'Y^H 1^ ((a + m ± + A' + {2n + ifr^^T')] . (4-4) 

± n=— oo 

If one takes only the lowest frequencies (ra = and n = —1) into account, Eq. fl4-4p becomes 

Q = 2G\a^ + Z\2) - 2f3-' In {{a + m ± fif + + n^T^) - 2/?-^ \n/3\ (4-5) 

In comparison to the KTV model Eq. fl2-14p . it is noted that the definition of temperature 
is different from that of the VJ model. Let us consider the high temperature limit, where 
the approximation in the KTV model with only the lowest Matsubara modes incorporated 
is expected to be justified. For this, the summation in Eq. (14-41) is carried out to yielcpS' 

n = 2G\a^ + Z\2) - 2T ^ log (e^^/^T ^ ^-e^/2T^^ ^4.5) 

± 

= 2G\a'' + A^)-Y, [E± + 2riog (l + e-^^±)] , (4-7) 
± 

where E± = a/(o" + ± /i)^ + A'^. By taking a limit T ^ 00 of Eq. (14 -6 1) and the effec- 
tive potential of the KTV model, Eq. (I2-I4p . we see that the VJ model is reduced to the 
KTV model at high temperatures, as expected, with the following identification of the two 
temperatures 

4Tvj = T|tv, (4-8) 

where Tyj (Tktv) denotes temperature in the VJ (KTV) model. We have checked this 
correspondence of temperature by looking at the values of chiral condensates a. At Tktv = 
2.0 (Tvj = 1.0), for example, a ratio of cr of the VJ model to that of the KTV model is 1.030 
at /i = 0.5 and 1.713 at /i = 1.0. At Tktv = 6.0 (Tyj = 9.0), they are 1.021 {fi = 0.5) and 
1.086 [fi = 1.0). As T increases, the difference between the VJ and the KTV models become 
smaller, and the range where the ratio is approximately 1.0 becomes wider near /i = 0. 

By numerically solving the saddle point equations for m = 0.1, a similar phase structure 
to that in the right panel of Fig. [1] is obtained in the real fi region, but with different 
temperature scale. The diquark condensate phase is located in the region 0.245 < /i < 0.519 
at Tvj = 0. The location, n = 0.245, of the lower critical point is in fairly good agreement 
with a lowest perturbative calculation with a small parameter m, fi^ = m/{2G'^){fic = 0.223). 
This region closes itself at Tyj = 0.194(= To). A cross over line starts at (/i,Tvj) = 
(0.34, 0.194) and moves towards higher T and smaller fi region as shown in Fig. [9], crossing 
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the T-axis at /i = and Tyj = 0.35(= Tco). It is noted that if one uses the relation Eq. (14- 8p . 
these values are in good agreement with the corresponding temperature of the KTV model; 
{Tof /A = 0.89174 = 0.198, and {Tcof /A = 1.18V4 = 0.348. In the imaginary fi = icp 




Fig. 9. Phase structure m fj, — T plane. In the imaginary fi region, additional chiral susceptibility 
peaks appear due to the periodicity 27r in 6* = (f)/T. m = 0.1. Solid and dotted lines indicate 
the second order critical line and cross over, respectively. Negative side of the horizontal axis 
indicates the imaginary chemical potential, 



Table VII. Values of data and the results of various fits at real values of /i for Tyj = 0.15 in the 
VJ model. The range of used data is 0.0 < < 0.3 and A(j) = 0.01. The lower (higher) critical 
point of the diquark condensation exists at = 0.270(0.472). 





fi = 0.0 


fi = 0.1 


H = 0.2 


H = 0.3 


fi = 0.4 


fi = 0.5 


data 


0.47942 


0.47414 


0.45345 


0.29623 


0.10337 


0.03243 


4th degree 


0.47955 


0.47473 


0.45812 


0.42327 


0.35944 


0.25159 


6th degree 


0.47944 


0.47426 


0.45512 


0.40878 


0.30490 


0.08376 


8th degree 


0.47943 


0.47417 


0.45409 


0.40047 


0.25671 


-0.13064 


10th degree 


0.47942 


0.47415 


0.45384 


0.39884 


0.24591 


-0.21725 


ratio 


0.47938 


0.47387 


0.45062 


0.34727 


0.99233 


0.63672 



region, the VJ model respects 27r symmetry in ^ = 0/T, in contrast to the KTV model. 
This is seen, when one sets Z\ = in the imaginary fi region, in the the effective potential 



2GV-2^1og exp 



a + m 9 

— 7^ ± 
2T 2 



+ exp 



0" + m 9 
=F ^- 

2r 2 



(4-9) 
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Table VIII. Values of data and the results of various fits at real values of fi for Tyj = 0.3 in the 
VJ model. The range of used data is 0.0 < (f) < 0.3 and A(j) = 0.01. A susceptibility peak is 
located at fj, = 0.245. 





= 


0.0 


/i = 0.1 


fi = 0.2 


fi = 0.3 


fi = 0.4 


fi = 0.5 


data 


0.2^ 


^040 


0.26792 


0.23090 


0.17619 


0.12217 


0.08200 


4th degree 


0.2^ 


^029 


0.26739 


0.22674 


0.15250 


0.0349 


-0.13954 


6th degree 


0.2^ 


^038 


0.26776 


0.22908 


0.16381 


0.0775 


-0.0853 


8th degree 


0.2^ 


^040 


0.26790 


0.23076 


0.17741 


0.1564 


0.34229 


10th degree 


0.2^ 


^040 


0.26792 


0.23113 


0.18232 


0.20021 


0.62457 


ratio 


0.2^ 


^026 


0.26725 


0.22569 


0.14681 


0.01016 


-0.23318 



due to the summation J2± '^^^^ positive and negative Matsubara frequencies, where T = Tyj. 
The cross over Une stated above crosses the T-axis at T^o and rises almost linearly. In 
addition to this, pairs of cross over lines appear, starting at (0 = 2mi, Tyj = Tco). In the low 
temperature region, the susceptibility peak appeared in the KTV model (subsction 12.2.21) . 
although it is very broad and low, does not appear. So, this structure of the peak is due to 
the approximation limiting the Matsubara frequencies. 

For the purpose of the analytic continuation, we studied the model by dividing the 
temperature into the three regions as in the previous section: Tyj < To (region I), To < 
Tyj < Tco (region II) and Tco < Tyj (region III). At temperatures in region I, the analytic 
continuation hits the lower critical point of the diquark condensate phase at /i = yUci- Up to 
this point, both of the polynomial and ratio types of fit work but beyond it the fits break 
down. The results for Tyj = 0.15 is listed in Table IVlTl At this temperature, /id = 0.270. 
At /X = 0.25, 0.27, 0.29, for example, the relative difference between the data and the 10th 
degree polynomial (ratio) type of fit are 0.58% (2.19%), 1.48% (3.54%), 23.15% (11.76%). 
In region II, we chose Tyj = 0.3. At this temperature, a susceptibility peak is located at 
/i = 0.245. In the similar way to the KTV case, both of the polynomial (with degree 4 to 10) 
and the ratio types of fit work well up to a point slightly higher than the cross over point, but 
beyond this point all extrapolations rapidly deviate from the data. The results for Tyj = 0.3 
is listed in Table IVIIII In contrast to these results, quite different results are obtained in 
region III. Figure [10] indicates the behaviors at Tyj = 0.5. See also Table IIX[ Although 
the polynomial type of fits show very slow convergence, the ratio type works very well, in 
agreement with the results in the KTV model in the previous section. The truncation of the 
Matsubara sum does not affect the properties of the chiral condensate at temperatures higher 
than the pseudo critical temperature. This fact may suggest that it is rather general feature 
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in accordance with the resuhs of MC simulations. ' The values of the fitting parameters 
used for all the temperatures are summarized in Table 1X1 

@ 

(m,Tvj)=(0.1,0.5) 

0.5 
0.4 

o 

0.3 
0.2 
0.1 





-0.4 -0.2 0.2 0.4 0.6 0.8 1 



Fig. 10. Analytical continuation of chiral condensate at Tyj = 0.5 in the VJ model for m=0.1. 

Table IX. Values of data and the results of various fits at real values of fi for Tyj = 0.5 in the VJ 
model. The range of used data is 0.0 < (p < 0.3 and A(p = 0.01. 





fi = 0.0 


fi = 0.2 


/i = 0.4 


fi = 0.6 


fi = 0.8 


fi = 1.0 


data 


0.09747 


0.09064 


0.07403 


0.05513 


0.03231 


0.02671 


4th degree 


0.09748 


0.09078 


0.07644 


0.07166 


0.10515 


0.21708 


6th degree 


0.09747 


0.09063 


0.07356 


0.04831 


-0.00952 


-0.19261 


8th degree 


0.09747 


0.09064 


0.07406 


0.05638 


0.05750 


0.17265 


10th degree 


0.09747 


0.09064 


0.07406 


0.05637 


0.05740 


0.17187 


ratio 


0.09747 


0.09061 


0.07374 


0.05391 


0.03589 


0.02131 



§5. Summary 

We studied analytic continuation in two-color QCD in terms of the chiral RMT with 
the lowest Matsubara frequencies as the temperature effect. In each of the three different 
temperature regions, we used polynomial and ratio types of fitting functions for the chiral 
condensate. In the temperature region higher than the pseudo critical point at /i = 0, the 
ratio type of fit works well, while the power types show slow convergence. This is the same 



data • 
power 4 — 

power 6 

power 8 

ratio 




'''.I 
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Table X. Parameters of various fits of the chiral condensate in the VJ model in accordance to 
a{(j)) = A + BcfP' + C(/>^ + L>(/>6 + E(t)^ + F^^o (polynomial; p) or a^cj)) = ^±§^ (ratio; r). Note 
that blank columns stand for terms not included in the fit. 



Tvj 


Fit 


A 


B 


C 


D 


E 


F 


d.o.f. 


xVd.o.f. 




P 


0.47955 


-0.46416 


-1.79061 








28 


1.23052 X 10-8 




P 


0.47944 


-0.49064 


-2.65467 


-6.8544 






27 


1.27281 X 10-1° 


0.15 


P 


0.47943 


-0.49517 


-2.92689 


-11.9795 


-29.7373 




26 


1.54657 X 10-12 




P 


0.47942 


-0.49592 


-2.99819 


-14.3132 


-60.5581 


-141.322 


25 


2.02449 X 10-1^ 




r 


0.47938 


-3.984 


-7.24552 








28 


2.44716 X 10-9 




P 


0.28029 


-1.27399 


-1.62144 








28 


7.7485 X 10-9 




P 


0.28038 


-1.25332 


-0.94699 


5.35106 






27 


3.38165 X 10-1° 


0.3 


P 


0.28040 


-1.24591 


-0.50153 


13.7362 


48.6624 




26 


1.36635 X 10-12 




P 


0.28040 


-1.24521 


-0.43433 


15.9355 


77.7099 


133.191 


25 


1.03238 X 10-1^ 




r 


0.28026 


-1.70342 


-1.50248 








28 


1.54665 X 10-8 




P 


0.09748 


-0.17932 


0.29892 








28 


3.40922 X 10-11 




P 


0.09747 


-0.18072 


0.25322 


-0.36259 






27 


1.35924 X 10-1^ 


0.5 


P 


0.09747 


-0.18067 


0.25605 


-0.30932 


0.30912 




26 


1.91676 X 10-20 




P 


0.09747 


-0.18067 


0.25605 


-0.30933 


0.30898 


-0.00063 


25 


1.99022 X 10-20 




r 


0.09747 


-0.04686 


1.37552 








28 


1.554 X 10-12 



feature as that found in the study of MC simulations.'^ In order to check whether this 
result may be affected by an approximation involved in this RMT model, we also studied 
a RMT model incorporating all the Matsubara frequencies^ and found that the same 
result is obtained in this temperature region. It might thus be rather general behavior, and 
reflect some non trivial dynamics of the quark gluon plasma with strong correlations.'^^' It is 
interesting that study in the imaginary region may also provide useful information on the 
phase structure in the real /i region.'^Sl' It would then be encouraging to study three-color 
QCD in terms of RMT. The results of this analysis will be reported in the forthcoming 
paper. We also investigated analytic continuation of the pseudo critical line. Its behavior is 
also similar to that in the MC simulations.'^ 

In view of the phase structure, some drawbacks of the RMT models should be mentioned. 
The first is a tricritical point of the diquark condensate transition,'l3'''2D which is missing 
in the RMT models. Because of their mean field nature, the first order phase transition at 
higher temperature side along the critical line is not detected. This concerns our region I, 
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where the analytic continuation breaks down at the critical line in any case, irrespective of 
the order of the transition. The second is a matter of the saturation of the interactions, which 
manifests as a decreasing behavior of the diquark condensation for /i > fid, as discussed in 
section [21 Temperatures in region III are outside of the region (region I) where the saturation 
concerns in our analysis, so we expect that the above stated behaviors in region III would 
be little affected, or that even if affected, it would occur in the much larger fi region. The 
last is the Roberge- Weiss symmetry, which reflects the center symmetry of the gauge group. 
This symmetry periodically induces critical lines at 6{= (p/T) = {k + l/2)7r(/c: integer). 
In the RMT model, this symmetry is missing. The existence of the critical line closest to 
the /i = axis affects the analytic continuation.'^^ In terms of the imaginary jj, variable, 
= T'k/2, the distance from the /i = axis is proportional to T, and its effect becomes less 
serious at higher temperatures. This may be a reason why the RMT model leads to the same 
conclusions as the MC simulations at higher temperature regions. In the lower temperature 
region, on the other hand, missing RW symmetry may cause a discrepancy between the 
lattice gauge theory and the corresponding RMT model. In fact, a use of a periodic fitting 
function is suggested in the low temperature region in the study of MC simulations.'^ This 
is, however, beyond the scope of the present analysis, unless some RMT model with the RW 
symmetry is considered. 
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